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Spin fluctuations in the stacked-triangular antiferromagnet YMnOj 
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The spectrum of spin fluctuations in the stacked-triangular antiferromagnet YMn03 was studied above 
the Neel temperature using both unpolarized and polarized inelastic neutron scattering. We find an in-plane 
and an out-of-plane excitation. The in-plane mode has two components just above Tjv, a resolution-limited 
central peak and a Debye-like contribution. The quasi-elastic fluctuations have a line-width that increases 
with q like Dq z and the dynamical exponent z=2.3. The out-of-plane fluctuations have a gap at the magnetic 
zone center and do not show any appreciable g-dependence at small wave-vectors. 

PACS: 75.40.Cx, 75.40.Gb, 78.70.Nx 



1. INTRODUCTION 

YM11O3 belongs to the family of RMn0 3 (R=Rare- 
carth) manganite ferroelectric compounds that crystal- 
lize in the hexagonal space-group P63C7T1 below the 
paraelectric-ferroelectric phase transition (~ 900 K). In 
YMnC>3 the Mn 3+ -ions form triangular layers well sep- 
arated from each other by Y-layers. Because in the fer- 
roelectric phase the lattice is distorted, the Mn-ions are 
slightly trimcrized. The large separation between adja- 
cent layers suggests that YMnC>3 forms a good candi- 
date of a geometrically frustrated 2-Dimensional anti- 
ferromagnet. 

The magnetic structure of hexagonal YMnC>3 was 
first investigated by Bertaut and Mercier and re- 
investigated later in more details by Munoz et al. [2]. 
Below Tjv ~70 K, the S=2 magnetic moment of Mn- 
ions are arranged in a 120° magnetic structure with the 
triangular layers at 2=0 and 2=1/2 being antiferromag- 
netically coupled. At saturation the magnetic moment 
is /i=2.9 fiB, i-e. significantly reduced from the expected 
4/ig of Mn 3+ spins, which was taken as evidence that 
even in the ordered phase strong spin fluctuations are 
present as a consequence of geometrical frustration [3]. 
Analysis of the spin- wave spectrum have confirmed the 
2-Dimensional character of the magnetic exchange inter- 
actions in YMnC>3 0] with the ratio of intra- to inter- 
plane exchange interactions being of the order ~2T0 2 . 
Whereas well-defined excitations are observed below the 
ordering temperature, a broad inelastic signal as well as 
short-range correlations between the Mn magnetic mo- 
ments within the triangular layers persist well above 
Tjv 0. 

In this work we investigate the q and temperature 
dependence of the spin excitations in YM11O3 close to 



the Neel temperature. We were motivated by the fact 
that the nature of the phase transition of frustrated 
magnets is still not completely understood and that the 
critical properties of stacked triangular antifcrromag- 
nets have received special attention since Kawamura 
proposed that the critical exponents in these systems 
form a new universality class. Experimental confirma- 
tion of the new class of (chiral) exponents was found 
in CsMnBr3 by unpolarized 7; and polarized neutron 
scattering & measurements as well as in Ho 0. Second 
an anomaly in the dielectric constant e with the electric 
field applied in the ab plane was found in YMnC>3 [TO] 
at Tjv- The nature of the coupling between electric and 
magnetic properties in hexagonal manganites is a sub- 
ject of intense debate ^5- Clearly it is required to char- 
acterise the behavior of the magnetic fluctuations in the 
vicinity of Tjv to understand the possible relationship 
with the magneto-dielectric effect in YMnOs |12j . 



2. EXPERIMENTAL 

Polycrystalline YMn03 was prepared by a solid 
state reaction. Starting materials of Y2O3 and MnC>2 
with 99.99% purity were mixed and grounded and then 
treated at temperature 1000-1200 C in air during at 
least 70h with several intermediate grindings. The 
phase purity of the compound was checked with con- 
ventional x-ray diffractometer (SIEMENS D500). The 
powder was hydrostatically pressed in the form of rods 
(8 mm in diameter and ~60 mm in length). The rods 
were subsequently sintered at 1300 C during 30h. The 
crystal growth was carried out using an Optical Floating 
Zone Furnace (FZ-T-10000-H-IV- VP-PC, Crystal Sys- 
tem Corp., Japan) with four 1000 W halogen lamps as 
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Fig 1: Temperature dependence of the intensity of (1 
0) magnetic Bragg peak. The line is a fit to the data 
with a power law. 



a heat source. The growth rate was 1.5 mm/h and both 
rods (feeding and seeding rod) were rotated at about 20 
rpm in opposite directions to ensure liquid homogeneity. 
A mixture of argon with 2% of oxygen at 5.5 bar was 
applied during growing. The crystal has the shape of a 
rod with 6 mm diameter and 2 cm height and a mosaic 
spread better than 1°. 

The measurements were performed at triple- axis 
TASP located on the cold source of the neutron spalla- 
tion source SINQ. The sample was mounted inside an 
He-flow cryostat with the crystallographic axis a* and 
c* in the scattering plane. The spectrometer was op- 
erated with the energy of the scattered neutrons kept 
fixed at k/=1.47A — 1 for measurements with unpolar- 
ized neutrons. 80' Soller collimators were installed in 
the incident beam and before the analyzer and the de- 
tector. With that configuration the energy resolution 
at zero energy transfer is 180 fieV. To reduce both the 
background and contamination by higher wavelength 
neutrons a cold Be-filter was installed in the scattered 
beam. The inelastic polarized neutron measurements 
were performed at k/=1.5lA _1 along (1 ± q, 0,0.1) 
and at different temperatures above Tpj. To perform 
longitudinal-polarization analysis remanent supermirror 
benders |KSj were inserted after the monochromator and 
before the analyser. The orientation of the polariza- 
tion was chosen perpendicular to the scattering plane. 
Because magnetic fluctuations with polarization factor 
parallel to the neutron spin occur in the non-spin flip 
channel, the non-spin-flip data (NSF) contains the in- 
plane-fluctuations and in the spin-flip channel (SF) only 
out-of-plane fluctuations are present, as will be shown 
below. 
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Fig. 2: a) Elastic scans along [1,0, 1] that show the in- 
crease of diffuse magnetic scattering in YMnC>3 when 
the temperature approaches Tat. b) Inelastic scan that 
show the presence of the central peak in YMnC>3. The 
insert shows the difference between scans measured at 
(1,0,-3.8) and (1,0,-3.5) that emphasizes the resolution- 
limited component. 



3. RESULTS AND DISCUSSION 

Figure ^ shows the intensity of the (10 0) mag- 
netic Bragg reflection as a function temperature which 
mirrors the square of the staggered magnetization. 
The transition temperature, as determined by tak- 
ing the derivative of the magnetization curve |14| . is 
Tat=72.1±0.05 K. The intensity of the Bragg peak fol- 
lows a power law I oc \(T/T N - 1)| 2/3 with /?=0.187(2). 
This value is lower than reported for the case of typical 
stacked-triangular antiferromagnet RbNiCi3 (/3=0.28), 
CsNiCl 3 09=0.28), CsMnBr 3 (0.2K /3 <0.25) and close 
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Fig. 3: a) Neutron diffuse intensity in YMn03. The 
background was measured at T=300 K and has been 
subtracted, b) Temperature dependence of the inverse 
of the correlation length. The line is a fit to the data 
with a power law. 



to the critical exponent obtained in VCI2 (/3=0.20) |SJ 
E|. Typical two-dimensional XY antiferromagnets 
(BaNi 2 (P0 4 ) 2 ), or ferromagnets (Rb 2 CrCl 4 , K 2 CuF 4 ) 
have a magnetization exponent that corresponds to the 
expected theoretical value /3=0.23 ^5]. We note that 
/3=0.19 in the three-dimensional triangular Ising anti- 
ferromagnet pni 

Figure [2^ shows elastic scans along the [1 I] di- 
rection as a function of temperature. Close to the Neel 
temperature, correlations between adjacent hexagonal 
planes give rise to broad magnetic scattering along c* 
that eventually appear to condense into the magnetic 
Bragg peak at Tjv- From Fig.|2h we see that this diffuse 
scattering is static on the time-scale of our experiment 
and corresponds to slow fluctuations of in-plane char- 
acter, as will be shown below. With increasing temper- 
ature these correlations disappear and only scattering 
across the [1 Z]-direction show a peak in the neutron 



cross-section. Hence, well above Tn correlations per- 
sist only between Mn moments located in the hexag- 
onal plane, as shown in Fig. |3| Here we describe the 
line-shape of the diffuse intensity by a Lorentzian pro- 
file convoluted with the resolution function of the spec- 
trometer: 
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where Qo is the position of the magnetic rod in recip- 
rocal space; || and _L denotes direction along and per- 
pendicular to the magnetic rod; Ku the inverse of the 
correlation length £ between Mn-spins in the hexagonal 
plane. Close to the Neel temperature, the temperature 
dependence of £ behaves like £(T) = 0.038(±0.005) • 
(T - T N f with v = 0.57(±0.06) as shown in Fig. EJd. 
We turn now to the behavior of the paramagnetic fluctu- 
ations in YMn03. Figure 01 shows typical energy scans 
performed at Q=(l-q,0,-3.8) and T=73.6 K. The in- 
clastic cross-section for an unpolarized neutron beam is 
given by 



d 2 <7 

dfldio 



^2(S a p - Q a Qf))S al3 (Q, uj), 



(2) 



a/3 



where 6 a p is the Kronecker symbol, a, [3 are the Carte- 
sian coordinates x,y,z, (Q, u/) denote the momentum 
and energy transfers from neutron to sample and Q = 
Q/|Q|. The first term in equation (0) is a selection rule 
that implies that only spin components perpendicular 
to the scattering vector contribute to the neutron scat- 
tering cross-section. Hence, for scattering vectors with 
large /-components like Q = (1 ± q, 0, —3.8), the inelas- 
tic spectrum contains essentially paramagnetic fluctua- 
tions with in-plane character (||). In addition, close to 
the Neel temperature, inelastic scans through the mag- 
netic rod show the resolution-limited central peak de- 
scribed by Eq. ^ as shown in Fig. Thus, to analyze 
the data shown in Fig. 0] we modeled the inelastic in- 
tensity /(Q, uS) in the following way 



.(oj) + S X c(Q)) 



R(Q,uj) + Bck (3) 



where Si„ c , = A5(u>) refers to the resolution- 
limited incoherent scattering that was measured 
at high temperature and S is a scale factor. The 
neutron scattering function Sp ar a. (Q> w)> which 
is related to the imaginary part of the dynami- 
cal susceptibility through 7r(g/is) 2 S']J Qra .(Q, oj) = 
F 2 (Q)(1 - exp(-^/fc B r))- 1 3xll(Q,c t ;), describes 
the line-shape of the paramagnetic scattering as a 
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Fig. 4: Constant q-scans in YMn0 3 at T=73.6 K and 
79.8 K that show the in-plane fluctuations, as explained 
in the text. The lines are the result of fit to the data 
using Eq. El 

function of momentum (Q) and energy (hui) transfer, 
respectively. F(Q) is the magnetic form factor of Mn. 
In Eq. |3 the symbol ® stands for the convolution with 
the spectrometer resolution function R(Q,u>) ^2] and 
Bck denotes the background level. We find that a 
Debye-like quasi-elastic line-shape for the imaginary 
part of the dynamical susceptibility 

3x"(Qo + q,")="X(Qo + q) 2 ^ V2 (4) 

reproduces the data adequately. x(q) is the static sus- 
ceptibility as in Eq. ^ taken relative to the antifer- 
romagnetic zone center (Qo=(l,0,0)) and T(q) is the 
damping of the paramagnetic fluctuations. Figure 0] 
shows the results of fits to the data at T=73.6 K from 
which we extract that the damping of the in-plane fluc- 
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Fig. 5: Constant q-scans in YMn0 3 at T=73.6 K 
that show the in-plane fluctuations in the non-spin flip 
(NSF) and the out-of-plane component in the spin-flip 
channel (SF), respectively. The lines are the result of 
fit to the data as described in the text. 

tuations T(q) evolves like r(0, T) + (U50±90)q z (meV) 
with r(0, 73.6 K)=0.23 ±0.07 (meV), and z=2.26±0.07. 
In contrast with data taken around (l±g,0,3.8) the in- 
elastic spectra around (l±g, 0,0.1) cannot be fitted with 
Eq. 0] This is an indication that the spectrum of para- 
magnetic fluctuations in YMJ1O3 consists of 2 modes, an 
in-plane (||) as well as an out-of-plane component (_L). 
To separate the ||- from the _L-fluctuations, it is nec- 
essary to use polarization analysis. A typical inelastic 
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Fig. 6: Temperature dependence of the damping of the 
paramagnetic fluctuations with in-plane polarization as 
measured with polarized neutrons. The line is a fit to 
the data as explained in the text. 

spectrum measured with polarization analysis is shown 
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Table 1: Parameter values for the damping (r(q) = 
"/q z (meV)) and the critical exponent (2) of the ||- 
fluctuations in YMn03. 
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pol. 
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2.53 


0.25 




78 


1028 


313 
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tic at the magnetic zone center and do not show any 
g-dependence for small wave- vectors. 
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in Fig. [3] that reveals an out-of-plane excitation in the 
spin-flip channel. In contrast with the ||-fluctuations, 
the ^-component is inelastic and is best described by a 
damped-harmonic oscillator function 



DHO l-exp(-aj/T)(co 2 -Q^) 2 +co 2 ^' w 

where £l q is the renormalized frequency and j q the 
damping of the excitation. Within the precision of the 
measurements we find that f2 g =3.9 ±0.3 (meV) and 
7 q =2 ±0.5 (meV) in the range of momentum values 
reached in the present experiment (q <0.1). 

The analysis of the data in the non-spin-flip chan- 
nel yields the same (^-dependence for the || -fluctuations 
than was obtained with the unpolarized set-up. The 
results are summarized in Tabic The mean value 
for the dynamical exponent of the || -fluctuations in 
YMnC>3 yields z ^2.3 and does not agree with the 
theoretical dynamical exponent z=l for the classical 
2D-triangular antiferromagnet |18| . Also the value ob- 
tained in YMn03 is quite different from the dynamical 
exponents z—1.5 expected for the 3D-Heisenberg anti- 
ferromagnet and measured e.g. in RbMnF3 ^U] and 
CsMnBr3 [201 ■ Finally we show the temperature depen- 
dence of r(0, T) in Fig.[|J]that increases almost linearly 
above Tat, r(0, T) = 0.4 + 0.07 * (T — T N ) 1A±0 - 2 . 

4. CONCLUSION 

To conclude, using both unpolarized and polarized 
inelastic neutron scattering, we showed that there are 
two magnetic excitations in the paramagnetic regime of 
YMnC>3 that have in-plane and an out-of-plane polariza- 
tion, respectively. The in-plane mode has a resolution- 
limited central peak and a quasi-elastic component with 
a line-width that increases as a function of momen- 
tum transfer q like oc q z with the dynamical exponent 
z=2.3. The presence of two time-scales in the spectrum 
of in-plane flucutations might be the signature for the 
coexistence of 2D and 3D fluctuations in the vicinity 
of T^r |221 ■ The out-of-plane fluctuations are inelas- 
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